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ABSTRACT 


An  experimental  program  was  carried  out  to  determine  the 
mass,  inertia,  and  stiffness  properties  of  a subsonic,  actively 
controlled,  variable  geometry,  wind  tunnel  model.  These  experimental 
measurements  were  then  employed  to  develop  a lumped  mass  dynamic 
model  from  which  vibrational  modes  and  frequencies  could  be  determined. 
The  analytical  modeshapes  and  frequencies  of  this  model  were  verified 
experimentally  for  one  of  the  geometry  configurations. 

An  analytical  analysis  was  then  conducted  for  the  uncontrolled 
model  to  predict  the  flutter  speeds  and  frequencies  during  wind  tunnel 
testing.  In  addition,  the  effect  on  flutter  speeds  and  frequencies 
of  various  mass  and  stiffness  changes  was  investigated  In  order  to 
specify  a final  test  configuration  which  would  remain  within  the 
model's  operating  limitations. 
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I . INTRODUCTION 


Active  control  technology  has  recently  become  an  integral 
part  of  the  design  process  for  advanced  aerospace  vehicles.  The 
application  of  this  new  technology  to  high  performance  aircraft  has 
necessitated  Investigation  of  new  methods  for  analyzing  and  experi- 
mentally testing  flutter  suppression  and  load  alleviation  systems. 

The  Increased  emphasis  of  control  innovations  being  employed  in 
Control  Configured  Vehicles  (CCV)  has  Indicated  that  dynamic  stability 
requirements  may  be  achieved  through  the  use  of  active  feedback  con- 
trol systems.1  Flutter  suppression  of  flexible  lifting  surfaces  Is  a 
natural  extension  of  the  new  CCV  concept  and  Initial  analytical 

O *>  A 

Investigations  * **  have  Indicated  that  feedback  control  laws  could 
be  developed  to  Increase  flutter  speeds  for  both  Isolated  lifting 
surfaces  and  aircraft  wings  with  external  stores.  More  recent 
experimental  verification  of  flutter  suppression  systems  has  been 
shown  In  Reference  5. 

With  the  support  of  the  Air  Force  Office  of  Scientific 

Research  and  the  assistance  of  the  Air  Force  Flight  Dynamics  Lab 

(AFFDL),  a graduate  level  research  program  was  begun  at  The  University 

of  Texas  to  study  active  flutter  suppression  on  Interfering  lifting 

surfaces  using  the  latest  and  most  accurate  unsteady  aerodynamic 

computer  programs.1,6,7,8  After  the  research  of  Griffin7  on  a small 

semirigid  model  with  three  degrees  of  freedom,  flutter  suppression 

1 8 

control  laws  were  developed  and  analyzed  by  Cwach  and  Plnnamanenl 
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for  the  more  complicated  variable-geometry  AFFDL  wing-tall  flutter 
model  shown  In  Fig.  1. 

A wind  tunnel  test  program  was  then  undertaken  to  experi- 
mentally verify  the  effectiveness  of  an  active  control  system  for 
this  model.  The  first  stage  of  the  test  program  consisted  of  the 
design  and  construction  of  the  electronic  and  hydraulic  control 
components  which  would  be  used  to  Implement  the  flutter  suppression 
feedback  control  laws.®  In  addition,  the  AFFDL  model  was  modified 
and  constructed  with  a remote  wing  sweep  mechanism  and  a hydraulically 
actuated  stabllator  control  surface  as  shown  In  Fig.  2. 

The  current  study  was  undertaken  to  Identify  the  system 
dynamic  parameters  and  to  Investigate  and  tune  the  passive  flutter 
properties  for  this  actively  controlled  model  prior  to  wind  tunnel 
testing.  The  tuning  procedures  utilized  mass  redistribution  and 
simple  stiffness  modifications  which  reduced  the  flutter  speeds  and 
frequencies.  The  system  dynamic  parameters  are  required  as  input  to 
the  active  control  law  determination  procedures.^  Additionally,  the 
passive  flutter  analysis  results  will  be  used  to  evaluate  model 
performance  during  the  wind  tunnel  testing  phase. 


II.  DYNAMIC  SYSTEM  IDENTIFICATION 


The  objective  In  defining  the  dynamic  system  was  to  use 
experimental  procedures  to  determine  an  accurate  mass  and  stiffness 
model  from  which  analytical  vibration  modes  and  frequencies  could 
be  calculated.  Experimental  verification  of  these  modes  and  frequen- 
cies would  then  be  attempted  for  the  60-degree  sweep  configuration. 
The  following  slmpllcatlons  In  the  analytical  modeling  were  employed 
to  facilitate  this  verification: 

(1)  The  effects  on  model  stiffness  related  to  pressure 
In  the  hydraulic  supply  lines  was  neglected. 

(2)  The  hydraulic  cylinders  were  assumed  to  be  rigid 
structural  links. 

These  simplifications  were  accounted  for  in  the  experimen- 
tal tests  by  disconnecting  the  hydraulic  supply  hoses  and  rigidly 
locking  the  cylinders  In  place. 

A.  Flutter  Model  Description 

The  actively  controlled  model  used  in  this  study  Is  a 
modified  version  of  a half-span,  variable-geometry,  passive  flutter 
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model  designed  by  the  Air  Force  Flight  Dynamics  Laboratory  and 
consists  of  a wing,  horizontal  tail,  and  fuselage  mechanism.  The 
original  model  was  modified  to  have  a remotely  operated  wing  sweep 
mechanism  and  a hydraulically  actuated  stabilator  control  surface 
(see  Fig.  1).  The  supporting  mechanism  of  Fig.  1 is  encased  in  an 
aerodynamic  housing  which  is  mounted  to  the  wind  tunnel  celling 
(see  Fig.  3).  Since  the  model  is  suspended  on  bearings,  it  is  free 
to  roll  under  the  influence  of  gravity.  That  is,  it  possesses  a 
rigid  roll  degree  of  freedom. 

The  fuselage  mechanism  and  wing  support  structure  was  designed 
to  allow  easy  variation  of  fuselage  stiffness,  wing-tail  horizontal 
separations,  and  wing  sweep  angle.  The  wing  spar  is  attached 
to  the  fuselage  roll  beam  by  a wing  carry-through  mechanism  which 
contains  both  roll  and  pitch  springs  (see  Fig.  2).  The  roll  spring 
Is  positioned  parallel  to  the  roll  axis  while  the  pitch  spring  is 
swept  20  degrees  from  vertical  (parallel  to  the  elastic  axis  for 
25-degree  wing  leading  edge  sweep).  The  tail  spar  Is  attached  to  the 
aft  fuselage  roll  beam  by  a rigid  carry-through  structure  which  is 
fitted  with  a bearing  to  allow  the  stabilator  to  pitch  about  a verti- 
cal axis.  Forward  and  aft  fuselage  sections  are  coupled  by  a torsion 
bar  whose  stiffness  is  varied  by  changing  the  spacing  of  the  fuselage 
torsion  clamps.  The  torsion  bar  Interconnector  (see  Fig.  2)  is 
not  Included  in  measurements  of  the  effective  spring  length. 


Geometry  modifications  Include  horizontal  wing-tall  separation 
(obtained  by  sliding  the  tall  support  along  the  aft  fuselage  roll  beam) 
and  variable  wing  sweep. 

The  wing  construction  consists  of  seven  wooden  aerodynamic 
sections  with  removable  panels  which  are  attached  to  the  tapered 
aluminum  spar  near  the  center  line  of  each  section  (see  Fig. 4 ). 

Hence,  the  elastic  axis  Is  coincident  with  the  spar  center  line.  The 
spar  is  located  at  37.5%  of  the  streamwise  chord  for  the  25-degree 
wing  sweep  configuration  and  at  an  angle  5 degrees  less  than  the 
leading  edge  sweep.  Each  section  is  separated  by  a slot  Into  which  is 
placed  a thin  piece  ' f *m  rubber  to  provide  aerodynamic  sealing. 

The  stabilator  is  of  identical  construction,  except  it  contains  only 
five  wooden  aerodynamic  sections.  The  tail  spar  is  located  at  the 
37.5%  stabilator  chord  line. 

B.  Experimental  Measurement  of  Mass,  Inertia,  and  Stiffness 

Before  assembling  the  wing  and  tall,  the  mass,  mass  static 
unbalance,  and  moment  of  inertia  about  the  spar  center  line  were 
determined  experimentally  for  each  wooden  aerodynamic  section.  From 
measurements  of  spar  dimensions  at  each  section  the  mass  and  Inertia 
values  for  the  spar  were  calculated  and  added  to  the  section  values. 

The  results  of  these  measurements  for  the  wing  and  tail  sections  and 
the  flexible  wing  carry-through  structure  are  shown  In  Table  I and 


TABLE  I 


EXPERIMENTAL  WING  MASS  DATA 


Section 

Location* 

Mass  x 102 

Unbalance  x 103 

Inertia  x 103 

(ft) 

(slug) 

(slug-ft) 

(slug-ft2) 

1 

.229 

2.390 

.782 

3.742 

2 

.625 

2.194 

.963 

3.268 

3 

.969 

1.881 

.671 

2.414 

4 

1.354 

1.556 

.604 

1.767 

5 

1.719 

1.299 

.432 

1.191 

6 

2.094 

1.103 

.257 

.789 

7 

2.484 

.992 

.183 

.541 

FI exlble 
Carry- 

8.011 

15.260 

8.157 

Through 

Section 

TABLE  II 

EXPERIMENTAL  TAIL  MASS  DATA 
Location*  Mass  x 102  Unbalance  x 103 

Inertia  x 103 

(ft) 

(slug) 

(slug-ft) 

(slug-ft2) 

1 

.198 

2.027 

1.123 

4.664 

2 

.552 

1.632 

1.019 

2.451 

3 

.917 

1.323 

.735 

1.480 

4 

1.271 

.880 

.394 

.597 

5 

1.620 

.524 

.083 

.216 

*D1stance  measured  along  the  spar  from  the  root. 


Table  II. 

The  mass  static  unbalance  for  each  section  about  the  spar 
center  line  was  determined  by  supporting  the  section  on  a knife-edge 
directly  beneath  the  spar  and  measuring  the  force  exerted  at  a known 
distance  from  the  support  (see  Fig.  5).  Inertia  measurements  for 
each  section  were  performed  by  supporting  the  section  as  above  and 
connecting  the  trailing  edge  to  two  springs  of  known  stiffness  as 
shown  In  Figure  6.  Next,  the  period  of  vibration  of  the 
section  was  measured  using  a Bentley  proximity  sensor  and  a frequency 
counter.  The  Inertia  could  then  be  calculated  using  the  following 
formula  which  Is  derived  from  the  equations  of  motion  for  a rigid 
body. 


°2(ki  + kj  Uy 

I = ! — - — - m D2  - Mx2 

S 2 2 s 

0)  0) 


(2.1) 


where 


Is  the  section  Inertia  about  the  spar  (lb-1n-sec2) 


Is  the  distance  from  the  support  point  to  the 
spring  attachment  (In) 

k^,  k2  are  the  support  spring  stiffnesses  (lb/in) 

Is  the  frequency  of  oscillation  (rad/sec) 

Is  the  weight  of  the  section  (lbs) 


■ VS  f. 


x Is  the  distance  of  the  center  of  gravity  above 

the  pivot  point  (In) 

ms  Is  the  mass  of  the  spring  attachment  device 

(lb-sec2/1n) 

M Is  the  mass  of  the  section  (1b-sec2/1n). 

The  Inertia  of  the  flexible  wing  carry-through  about  the  pitch  spring 
axis  was  measured  In  the  same  manner. 

Moments  of  Inertia  about  the  model  roll  axis  were  also 
measured  for  the  forward  and  aft  fuselage  sections  and  are  listed  In 
Table  III.  These  measurements  were  accomplished  by  accurately 
measuring  the  pendulum  frequency  of  the  component  and  then  measuring 
the  angular  deflection,  0,  for  a known  applied  force,  F,  as  shown  In 
Fig.  7.  Provided  that  the  deflection  angle  Is  kept  small  and  the 
mirror  Is  positioned  close  to  the  pivot  axis,  the  roll  Inertia  can 
be  calculated  by  the  following  formulas. 


where 

dg  Is  the  laser  beam  deflection  (In) 

r Is  the  distance  over  which  the  laser  beam  Is 

deflected  (ref.  Fig.  7) 


< iTTe  > 1 


(2.3) 


where 


Is  the  applied  force  (lbs) 


Is  the  pendulum  frequency  (rad/sec) 


Is  the  deflection  angle 


Is  the  distance  from  the  roll  axis  to  the 
force  application  point  (Ft) 


With  the  model  completely  assembled,  the  roll  frequencies 
were  measured  for  wing  sweeps  of  25,  45,  and  60  degrees  and  are 
presented  In  Table  IV.  All  of  the  above  measurements  were  made  with 
the  hydraulic  supply  hoses  disconnected. 

Elastic  properties  for  the  wing  and  tall  spars  were  deter- 
mined by  measuring  the  bending  and  torsional  Influence  coefficients 
with  the  root  effectively  cantilevered.  The  bending  coefficients 
were  measured  directly  with  an  LVDT  that  gave  a voltage  readout  on  an 
accurate  Integrating  digital  voltmeter  (see  Figures  8 and  10).  The 
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TABLE  III 

FUSELAGE  ROLL  INERTIA 


Component 


Ie  x 10  (slug-ft2) 


♦Forward  Fuselage  .623 

♦♦Aft  Fuselage  .275 


♦Forward  Fuselage  - Includes  leading  edge  glove,  wing  carry-through 

structure  and  all  components  of  the  fuselage 
mechanism  forward  of  the  torsion  spring 
Interconnector. 

♦♦Aft  Fuselage  - Includes  tall  carry-through  structure  and  all 

components  of  the  fuselage  mechanism  aft  of  the 
torsion  spring  Interconnector. 


TABLE  IV 

MODEL  PENDULUM  (ROLL)  FREQUENCIES 


Wing  Sweep 


fR  (Hz) 


25° 

45° 

60° 


0.78 

0.84 

0.87 


TABLE  V 

WING  CARRY-THROUGH  STRUCTURE  INFLUENCE  COEFFICIENTS 
Croll  * 5,8518  x 10’%  (rad/ln-lb) 

Cp1 tch  " 4*9806  x 10~ " ( rad/ in-1 b) 
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deflection  was  taken  at  each  of  the  spanwlse  stations  for  loads 
applied  at  all  the  other  stations.  The  voltage  was  converted  directly 
to  deflections  by  using  the  LVDT  calibration  factor.  Torsional 
Influence  coefficients  were  obtained  by  measuring  the  rotational 
deflection  at  a station  due  to  an  applied  moment  at  that  point  (see 
Figures  9 and  11).  An  optical  leverage  technique10  gave  deflections 
of  a laser  beam  which  could  be  converted  directly  to  rotational 
deflections  with  Eq.  (2.2).  All  bending  and  torsional  Influence 
coefficients  were  measured  for  several  applied  loads  and  moments 
and  then  averaged  over  all  loads  and  reciprocal  deflection  loading 
stations.  The  wing  and  tall  bending  and  torsional  Influence  coeffi- 
cients are  presented  In  Tables  VI  and  VII. 

Pitch  and  roll  Influence  coefficients  of  the  wing  carry- 
through  structure  were  determined  by  applying  a moment  to  each  spring 
and  measuring  the  rotational  deflection  of  the  spring  with  an  optical 
leverage  technique.  These  data  are  given  In  Taole  V. 

Definitions  of  fuselage  torsional  stiffness  versus  torsion 
bar  length  were  obtained  with  the  fuselage  mechanism  assembled.  The 
aft  roll  beam  was  clamped  rigidly  and  load  Increments  were  applied 
to  the  forward  roll  beam  creating  a moment  on  the  torsion  bar.  The 
torsion  bar  length  was  changed  by  moving  the  torsion  bar  clamps.  A 
mirror  was  attached  to  the  forward  roll  beam  and  the  deflections  of 
a laser  light  source  were  recorded  for  several  load  Increments  and 
converted  to  rotational  Increments  with  Eq.  (2.2).  Torsional  stiffness 
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TABLE  VI 

WING  INFLUENCE  COEFFICIENTS 


Wing  Bending  Along  Spar  (In/lb) 


.0004 

.0011 

.0018 

.0025 

.0032 

.0039 

.0047 

.0011 

.0056 

.0100 

.0149 

.0197 

.0244 

.0295 

.0018 

.0100 

.0214 

.0340 

.0466 

.0591 

.0722 

.0025 

.0149 

.0340 

.0615 

.0883 

.1158 

.1443 

.0032 

.0197 

.0466 

.0883 

.1385 

.1884 

.2409 

.0039 

.0244 

.0591 

.1158 

.1884 

.2748 

.3649 

.0047 

.0295 

.0722 

.1443 

.2409 

.3649 

.5176 

Wing  Torsion  About  Soar 

(rad/ In 

.00008 


.00022- 


.00046- 


00078- 


.00124- 


.00197- 


.00315 


TABLE  VII 


HORIZONTAL  TAIL  INFLUENCE  COEFFICIENTS 


Tall  Bending  Along  Spar  (In/lb) 


0004 

.0010 

.0015 

.0022 

.0027 

0010 

.0045 

.0085 

.0125 

.0165 

0015 

.0085 

.0199 

.0320 

.0444 

0022 

.0125 

.0320 

.0589 

.0872 

0027 

.0165 

.0444 

.0872 

.1446 

.00008- 


Tall  Torsion  About  Spar  (rad/ln-lb) 


.00026 


♦ 


00053 

nnim  

1 

• v/in  u i 

.00188 

■■■■■■■■■I 


•WsitSsSjw. 


mi 


was  then  computed  from 

✓ • 

V 


K = £ 

* A6 


(2.4) 


where 


K 

AT 


Is  the  torsional  stiffness  ( ‘p~jb) 


Is  the  Increment  of  moment  due  to  an  added 
weight  (In-lb) 


Ad 


Is  the  resulting  rotational  deflection  (rad). 


This  stiffness  Information  was  then  used  to  calculate  a 
curve  of  torsional  frequency  versus  torsion  bar  length.  The  uncoupled 
fuselage  torsional  frequency,  w0,  was  computed  by  assuming  that  the 
entire  wing  and  tall  units  were  rigid  and  acted  as  two  lumped  Inertias. 
The  vibrational  frequency  of  such  a two  degree  of  freedom  system 
which  has  a rigid  body  roll  mode  can  be  computed  from  (see  Ref.  11) 


where 


K 

Ir 


. . 7*515 

V 'f'a 


(2.5) 


ft-lb. 


Is  the  fuselage  torsional  stiffness  ( ) 

Is  the  Inertia  of  the  forward  fuselage  section 
and  wing  about  the  roll  axis  (lbs-sec2-ft) 
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is  the  inertia  of  the  aft  fuselage  section  and 
tall  about  the  roll  axis  (lbs-sec2-ft). 

The  curve  of  fuselage  torsional  frequency,  u)Q,  versus  torsion  bar 
length  for  a wing  sweep  of  60  degrees  1$  shown  In  Figure  16. 

C.  Analytical  Vibration  Analysis 

The  dynamic  model  from  which  the  natural  modes  and  frequencies 
were  calculated  was  constructed  from  the  experimental  measurements  of 
mass.  Inertia,  and  stiffness  described  In  the  preceding  section.  It 
consists  of  28  lumped  mass  points  and  a flexible  wing,  tall,  and  fuse- 
lage. As  outlined  In  Appendix  A,  a computer  code  determines  the  appro- 
priate mass  point  locations  and  assembles  the  model  flexibility  matrix 
for  the  required  wing  sweep  angle.  The  mass  modeling  for  the  60-degree 
sweep  configuration  Is  shown  In  Figure  12. 

Each  wing  sweep  configuration  In  the  vibration  analysis  Is 
characterized  by  Its  roll,  pitch,  and  fuselage  flexibilities,  and  by 
the  amount  of  lumped  mass  added  to  nodes  6 and  7 on  the  main  wing.  The 
ratio  of  the  coupled  wing  cantilever  bending  frequency,  w^,  to  the 

uncoupled  fuselage  torsional  frequency,  u)Q,  Is  used  as  a measure  of  the 

g 

coupling  between  wing  and  tall  and  can  be  computed  for  each  combination 
of  stiffness,  mass  and  wing  sweep  angle.  The  wing  cantilever  bending 
frequency  Is  the  first  flexible  mode  frequency  obtained  with  the  fuselage 
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rigid,  while  the  fuselage  torsional  frequency  is  determined  as  shown  in 
Section  II.B. 

Typical  results  for  the  computed  natural  frequencies  of  the 
first  six  modes  are  listed  In  Table  VIII  and  the  modeshapes  for  cases 
2 and  9 are  tabulated  In  Tables  IX  and  X.  Also,  three  dimensional 
computer  plots  of  the  modeshapes  for  several  of  the  configurations  are 
shown  In  Appendix  B.  The  configurations  listed  In  Table  VIII  were 
determined  by  modifications  made  to  the  dynamic  model  In  an  attempt 
to  reduce  flutter  speeds  and  frequencies.  An  explanation  of  these 
attempts  to  reduce  the  flutter  speed  will  be  given  in  the  section  on 
flutter  analysis. 

For  all  configurations  examined,  the  generalized  masses  and 
stiffnesses  were  calculated  from  the  computed  modes  and  frequencies  as 
shown  In  Appendix  A.  In  addition.  Inertia  terms  for  the  stabllator 
control  surface  were  defined  by^ 


ff  h (x,y)  h,  (x,y)  m (x,y)  dx  dy 
r 6j 


(2.6) 


or,  for  a lumped  mass  model, 

v • <w>]<v, 

rj  r j 

where 

(h«}T  Is  the  transpose  of  the  vector  of  deflections 
* r 

of  the  r*h  absolute  mode  (Ref.  Appendix  A) 
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{h^  Is  the  vector  of  deflections  on  the  j**1  control 

surface  due  to  a unit  deflection  of  the  jth  control. 
The  deflections  are  defined  over  the  area  of  the  jth 
control  surface  and  are  zero  on  the  remainder  of  the 
lifting  surface.  (For  the  stabllator  control  this 
is  the  vector  of  distances  from  the  pitch  axis 
for  a unit  rotation.) 

[M]  Is  the  lumped  mass  matrix. 

Typical  generalized  mass,  generalized  stiffness,  and  control  Inertia 
terms  for  configurations  2 and  9 of  Table  VIII  are  given  in  Table  XI. 

D.  Experimental  Vibration  Analysis 

An  evaluation  of  the  mathematical  modeling  was  carried  out  by 
experimentally  measuring  the  flexible  modes  and  frequencies  for  configu- 
ration 2 of  Table  VIII.  This  configuration  was  chosen  for  verification 
since  It  would  be  analyzed  thoroughly  in  Initial  wind  tunnel  tests  of 
the  model's  active  control  system.  Also,  only  the  first  four  flexible 
modeshapes  were  measured  although  frequencies  were  verified  for  five 
modes.  The  fifth  flexible  modeshape  was  not  measured  because  of  the 
very  small  deflection  values  at  this  frequency. 

The  experimental  apparatus  for  measuring  the  natural  modes  and 
frequencies  Is  shown  schematically  In  Figure  13.  In  this  arrangmment  the 
model  Is  excited  by  a shaker  attached  to  the  rigid  tall  carry-through 
structure.  The  model  Is  tuned  to  each  natural  frequency  by  displaying 


mm 
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TABLE  VIII 

NATURAL  VIBRATION  FREQUENCIES  FROM  EXPERIMENTAL 
MASS,  INERTIA  AND  STIFFNESS** 


No. 

Sweep 

n 

we 

Modification 

tMode  1 
(Hz) 

Mode  2 
(Hz) 

Mode  3 Mode  4 Mode  5 Mode  6 
(Hz)  (Hz)  (Hz)  (Hz) 

1 

O 

O 

AO 

.8 

None 

.87 

3.7 

8.1 

14.2 

19.0 

21.5 

2 

60° 

.6 

None 

.87 

4.0 

8.8 

14.3 

19.2 

23.2 

3 

60° 

.4 

None 

.87 

4.3 

9.8 

14.4 

19.2 

26.5 

4 

60° 

* 

Ml, (a.) 

.85 

3.7 

8.2 

14.1 

18.6 

22.3 

5 

60° 

* 

M2, (a.) 

.84 

3.5 

7.8 

13.9 

18.2 

21.9 

6 

O 

O 

AO 

* 

R»(b. ) 

.87 

3.8 

8.6 

14.1 

18.4 

22.9 

7 

60° 

* 

P,(c.) 

.87 

3.9 

8.5 

13.4 

19.1 

23.1 

8 

O 

O 

AO 

.6 

PM,(d. ) 

.84 

3.2 

7.0 

12.9 

17.8 

21.0 

9 

45° 

* 

None 

.84 

4.1 

9.7 

13.9 

19.2 

23.3 

10 

45° 

# 

PM,(d.) 

.80 

3.4 

7.7 

12.3 

18.0 

21.2 

11 

25° 

* 

None 

.78 

4.1 

10.4 

13.8 

19.1 

23.4 

12 

25° 

# 

PM,(d. ) 

.73 

3.5 

8.5 

11.8 

18.3 

21.3 

**  (see  also.  Appendix  B) 

+ Pendulum  frequencies  were  experimentally  measured  for  the  unmodified 
model  and  then  corrected  when  lumped  masses  were  added. 

* Same  fuselage  stiffness  as  configuration  No.  2 

# Same  fuselage  stiffness  as  configuration  No.  8 

(a.)  1.  10*  of  wing  mass  added  at  tip  trailing  edge 

2.  20*  of  wing  mass  added  at  tip  trailing  edge 

(b.)  Roll  spring  flexibility  Increased  by  50* 

(c.)  Pitch  spring  flexibility  Increased  by  50* 

(d.)  50*  Increase  In  pitch  spring  flexibility  and  20*  of  wing  mass 
added  at  tip  trailing  edge. 
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TABLE  IX 

ANALYTICAL  MODES  OF  MODIFIED  AFFDL  MODEL 
60°  Sweep  , — = 0.6 


Structural 

X 

Y 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Mode  6 

Control 

Point 

(FT) 

(FT) 

.87  Hz 

4.0  Hz 

8.8  Hz 

14.3  Hz 

19.2  Hz 

23.2  Hz 

1 

.41 

-.19 

-.073 

.013 

.022 

.213 

-.017 

-.055 

2 

.74 

.04 

-.103 

.046 

.013 

.317 

.043 

.007 

3 

1.00 

.26 

-.131 

.084 

-.012 

.375 

.090 

.070 

4 

1.30 

.50 

-.162 

.139 

-.080 

.387 

.106 

.109 

5 

1.58 

.74 

-.193 

.201 

-.176 

.307 

.062 

.056 

6 

1.87 

.98 

-.225 

.272 

-.303 

.140 

-.062 

-.118 

7 

2.17 

1.23 

-.257 

,351 

-.460 

-.094 

-.262 

-.386 

8 

-.00 

• 40 

-.150 

-.010 

.293 

-.087 

.007 

-.086 

9 

.33 

.62 

-.177 

• 023 

.294 

-.046 

.079 

.024 

10 

.62 

• 80 

-.201 

.063 

.267 

-.041 

.141 

.150 

11 

.95 

1.00 

-.227 

.118 

.197 

-.072 

.175 

.249 

12 

1.27 

1.19 

-.252 

.182 

.086 

-.161 

.144 

.256 

13 

1.59 

1.39 

-.277 

.255 

-.060 

-.317 

.027 

.136 

14 

1.92 

1.59 

-.304 

.335 

-.239 

-.521 

-.171 

-.106 

15 

.01 

-.14 

-.079 

-.016 

.118 

.017 

-.052 

-.101 

16 

-.46 

.05 

-.103 

-.011 

.170 

.008 

-.031 

-.099 

17 

-.17 

-.43 

-.042 

-.023 

.041 

.031 

-.085 

-.103 

18 

2.88 

-.15 

-.078 

-.082 

-.038 

.002 

-.082 

.097 

19 

3.05 

.18 

-.120 

-.130 

-.074 

-.008 

-.044 

.079 

20 

3.23 

.49 

-.162 

-.182 

-.123 

-.029 

.086 

-.030 

21 

3.39 

.82 

-.204 

-.237 

-.183 

-.059 

.296 

-.223 

22 

3.56 

1.14 

-.245 

-.293 

-.250 

-.096 

.569 

-.488 

23 

2.21 

.40 

-.149 

-.155 

-.069 

.008 

-.195 

.237 

24 

2.50 

.62 

-.177 

-.189 

-.098 

-.002 

-.154 

.223 

25 

2.77 

.87 

-.211 

-.232 

-.143 

-.023 

-.026 

.123 

26 

3.03 

1.11 

-.242 

-.276 

-.198 

-.054 

.197 

-.086 

27 

3.28 

1.36 

-.274 

-.323 

-.262 

-.091 

.488 

-.374 

28 

2.39 

0.00 

-.097 

-.100 

-.043 

.006 

-.133 

.153 
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TABLE  X 

ANALYTICAL  MODES  OF  MODIFIED  AFFDL  MODEL 
45°  Sweep + 


Structural 

Control 

X 

Y 

Mode  1 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Mode  6 

Point 

(FT) 

(FT) 

.84  Hz 

4.1  Hz 

9.7  Hz 

13.9  Hz 

19.2  Hz 

23.3  Hz 

1 

.44 

-.07 

-.079 

.005 

.050 

.235 

-.016 

-.079 

2 

.70 

.23 

-.114 

.033 

• 061 

.319 

.050 

.002 

3 

.90 

.51 

-.147 

.068 

• 042 

.363 

.105 

• 102 

4 

1.13 

.82 

-.183 

.117 

-.030 

.374 

.129 

.167 

5 

1.34 

1.12 

-.219 

.173 

-.145 

.312 

.084 

.107 

6 

1.55 

1.43 

-.255 

.238 

-.306 

.183 

-.042 

-.080 

7 

1.78 

1.75 

-.293 

• 309 

-.507 

.003 

-.233 

-.333 

8 

-.11 

.39 

-.133 

.000 

.271 

-.117 

.000 

-.084 

9 

.16 

.68 

-.167 

.028 

.292 

-.099 

.070 

.027 

10 

.39 

.93 

-.197 

.063 

.275 

-.107 

.128 

• 150 

11 

.66 

1.21 

-.230 

.111 

• 204 

-.137 

.156 

.238 

12 

.91 

1.48 

-.260 

.168 

.079 

-.206 

.117 

.216 

13 

1.17 

1.75 

-.292 

.233 

-.096 

-.321 

-.002 

.070 

14 

1.44 

2.03 

-.325 

• 305 

-.314 

-.467 

-.188 

-.166 

15 

.01 

-.14 

-.071 

-.014 

.111 

.021 

-.053 

-.109 

16 

-.46 

.05 

-.093 

-.007 

.155 

.017 

-.032 

-.105 

17 

-.17 

-.43 

-.038 

-.024 

.046 

.028 

-.084 

-.113 

18 

2.88 

-.15 

-.070 

-.089 

-.028 

.001 

-.083 

.100 

19 

3.05 

.18 

-.108 

-.143 

-.058 

-.009 

-.044 

.083 

20 

3.23 

.49 

-.145 

-.199 

-.103 

-.027 

.087 

-.026 

21 

3.39 

.82 

-.183 

-.260 

-.159 

-.054 

.299 

-.226 

22 

3.56 

1.14 

-.221 

-.322 

-.221 

-.087 

.573 

-.503 

23 

2.21 

.40 

-.134 

-.169 

-.049 

.005 

-.198 

.241 

24 

2.50 

.62 

-.160 

-.207 

-.075 

-.005 

-.156 

.231 

25 

2.77 

.87 

-.190 

-.254 

-.116 

-.023 

-.026 

.132 

26 

3.03 

1.11 

-.218 

-.303 

-.168 

-.050 

.201 

-.080 

27 

3.28 

1.36 

-.247 

-.355 

-.229 

-.083 

.494 

-.380 

28 

2.39 

0.00 

-.088 

-.110 

-.030 

.004 

-.135 

.156 

+ Same  fuselage  stiffness  as  60°  sweep  case. 
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TABLE  XI 

GENERALIZED  MASS  AND  STIFFNESS  MATRICES  OF  MODIFIED 
60°  Wing  Sweep  , -7-  = 0.6 


mn  = 8.678  x 10-3 
n»22  ■ 6.844  x 10-3 
m33  = 9.536  x 10'3 
m44  = 8.724  x 10-3 
nigg  * 7.558  x 10-3 
m66  = 1.006  x IQ-2 


= 8.436  x 10-3 
n>22  - 6.863  x 10"3 
m33  » 9.278  x 10"3 
m44  * 8.980  x 10-3 
m55  « 7.545  x 10"3 
m66  * 1.068  x 10'2 

* Stabllator  control  surface 
t Same  fuselage  stiffness  as 


kn  = .2593 

* 

m61 

k22  - 4.421 

m62 

k33  = 29.25 

m53 

k44  « 70.13 

m64 

k55  - 109.5 

m65 

k66  - 212.8 

m66 

45°  Wing  Sweept 

* 

kn  - .2350 

m61 

k22  “ 4.586 

m<52 

k33  - 34.14 

m63 

k44  = 68.10 

m64 

k55  - 109.5 

m65 

k66  - 229.7 

m66 

60°  sweep  case 


AFFDL  MODEL 


« -3.105  x 10-3 
= -3.524  x 10-3 
- -2.477  x 10-3 
= -6.496  x 10-4 
= 2.405  x 10'3 
* -1.272  x 10'3 


* -2.796  x 10‘3 

* -3.863  x 10'3 
■ -2.089  x 10'3 

* -6.084  x 10‘4 
= 2.436  x 10"3 

* -1.252  x 10'3 


All  Dimensions  In  Ft.  and  Lbs 
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0 

on  an  oscilloscope  the  displacement  sensor  output  versus  the  shaker  load 
voltage  to  form  a llssajou  figure.  With  this  procedure  the  resonant 
frequency  for  a mode  can  be  accurately  tuned  by  finding  the  point  of  90 
degree  phase  shift  between  the  forcing  function  and  the  structural 
displacement. 

The  modeshapes  were  measured  with  an  adjustable  displacement 

sensor  as  shown  In  Figure  14.  While  these  displacements  were  being 

measured,  an  additional  proximity  sensor  at  the  shaker  attachment  point 

was  monitored  to  maintain  a constant  excitation  amplitude.  This  proximity 

sensor  could  also  be  displayed  against  the  moveable  sensor  to  check 

relative  phasing  between  the  excitation  point  and  any  other  structural 
l 

location. 

As  discussed  In  the  Introduction  to  Section  II,  the  experi- 
mental modes  and  frequencies  were  Initially  measured  with  the  hydraulic 
supply  lines  disconnected  and  the  actuating  cylinders  locked  rigidly  In 
place.  The  experimental  modeshapes  were  then  normalized  (see  Appendix  A) 
and  are  listed  In  Table  XII  with  their  resonant  frequencies.  Appendix  C 
displays  computer  plots  of  the  measured  flexible  modes,  all  of  which  are 
In  very  good  agreement  with  the  computed  values.  The  calculated  frequen- 
cies are  all  within  5%  of  the  measured  values.  Indicating  that  the  elastic 
and  mass  properties  of  the  model  were  simulated  quite  closely. 


Natural  frequencies  for  all  six  modes  were  then  remeasured 
with  the  hydraulic  lines  connected  while  the  actuator  cylinders  remained 
locked.  These  frequencies  are  given  In  Table  XIII  and  are  slightly 
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TABLE  XII 

EXPERIMENTAL  MODES  OF  MODIFIED  AFFDL  MODEL 

wh 

60°  Sweep  , = 0.6 


Structural 

Control 

X 

Y 

Mode  2 

Mode  3 

Mode  4 

Mode  5 

Point 

(FT) 

(FT) 

4.2  Hz 

9.0  Hz 

15.1  Hz 

19.2  Hz 

1 

.41 

-.19 

-.015 

.014 

.122 

-.006 

2 

.74 

.04 

-.044 

.010 

• 188 

.041 

3 

1.00 

.26 

-.071 

-.005 

.233 

.031 

4 

1.30 

.50 

-.129 

-.046 

.229 

.031 

5 

1.58 

.74 

-.197 

-.120 

• 146 

.044 

6 

1.87 

.98 

-.293 

-.235 

.06? 

-.085 

7 

2.17 

1.23 

-.382 

-.360 

-.201 

-.150 

8 

-.00 

.40 

.025 

.257 

-.083 

• 006 

9 

.33 

.62 

-.014 

.269 

-.087 

.088 

10 

.62 

.80 

-.048 

.235 

-.097 

.110 

11 

.95 

1.00 

-.111 

.187 

-.097 

.13? 

12 

1.27 

1.19 

-.172 

.091 

-.059 

.119 

13 

1.59 

1.39 

-.229 

• 019 

-.135 

.081 

14 

1.92 

1.59 

-.322 

-.178 

-.302 

-.038 

15 

.01 

-.14 

• 018 

.074 

.021 

-.056 

16 

••46 

.05 

.036 

.149 

-.049 

-.056 

17 

-.17 

-.43 

.025 

.043 

.028 

-.078 

18 

2.88 

-.15 

• 125 

-.072 

.056 

-•?13 

19 

3.05 

• 18 

.122 

-.096 

-.010 

-.138 

20 

3.23 

.49 

.172 

-.168 

-.104 

-.025 

21 

3.39 

• 82 

• 236 

-.252 

-.264 

.229 

22 

3.56 

1.14 

.315 

-.351 

-.469 

.548 

23 

2.21 

.40 

• 136 

-.077 

• 111 

-.279 

24 

2.50 

.62 

.189 

-.130 

• 04? 

-.207 

25 

2.77 

.87 

• 229 

-.202 

-.094 

-.038 

26 

3.03 

1.11 

.279 

-.274 

-.274 

.229 

27 

3.28 

1.36 

.307 

-.360 

-.472 

.508 

28 

2.39 

0.00 

• 096 

-.050 

.059 

-.207 

I 

1 


1 
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TABLE  XIII 

EXPERIMENTAL  VIBRATION  FREQUENCIES  WITH 
HYDRAULICS  CONNECTED 

“h 

60°  Sweep  , ■ 0.6 


Mode  Frequency  (Hz) 


1 

2 

3 

4 

5 

6 


1.45 

4.2 

9.1 

15.8 

19.0 

22.0 


TABLE  XIV 

EXPERIMENTAL  STRUCTURAL  DAMPING  WITH 
HYDRAULICS  CONNECTED 


60°  Sweep 


% 

“6 


0.6 


Mode  Structural  Damping 


1 0.058 

2 0.013 

3 0.041 
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shifted  from  those  appearing  In  Table  XII,  Indicating  that  the  hydraulic 
lines  caused  an  Increase  In  stiffness  for  the  first,  third  and  fourth 
inodes  while  adding  an  apparent  mass  to  modes  five  and  six.  These  new 
modeshapes  appeared  to  be  similar  to  the  old  ones,  but  additional 
measurements  will  be  necessary  to  quantify  them  exactly. 

Damping  was  also  measured  for  the  first  three  modes  with  the 
hydraulic  lines  connected.  With  the  shaker  disconnected,  modes  one  and 
two  were  excited  by  hand.  The  method  of  logarithmic  decrement^’^  was 
then  used  to  obtain  damping  estimates  from  decay  traces  on  an  oscillo- 
scope. For  harmonic  motion  at  the  resonant  frequency 

* * 2T  ■ »r  • (sr)  <2-8) 

where 

g Is  the  structural  damping 

N Is  the  nuafcer  of  cycles 

XP  Is  the  Initial  amplitude 

XQ  Is  the  amplitude  after  n cycles. 

Damping  for  the  third  mode  was  obtained  by  exciting  the  model  with  the 
shaker  and  accurately  measuring  the  phase  lag  angle  of  the  response  for 
slight  frequency  variations  near  the  resonant  point.  The  variable 
phase  oscillator  shown  In  Figure  13  was  used  to  measure  the  phase 
variations.  Damping  values  could  then  be  calculated  from 
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where 
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g ■ tan  <J>  (1  - r2) 


(2.9) 


g Is  the  structural  damping 

♦ Is  the  phase  angle  by  which  the  response  lags  the 
forcing  function 

r Is  the  ratio  of  excitation  frequency  to  resonant 
frequency. 


The  results  of  the  damping  measurements  are  given  In  Table  XIV. 


III.  ANALYTICAL  FLUTTER  ANALYSIS 


The  objective  of  the  current  flutter  analyses  was  to  evaluate 
the  uncontrolled  flutter  speeds  of  the  modified  University  of  Texas 
version  of  the  AFFDL  model  and.  If  necessary,  reduce  these  speeds  and 
frequencies  to  a level  more  compatible  with  the  operating  limitations 
during  the  wind  tunnel  tests.  The  60  degree  sweep  configuration  was 
examined  for  various  combinations  of  fuselage  stiffness  and  wing-tail 
horizontal  separation,  X,  as  shown  In  configurations  1 - 6 of  Table  XV. 
Configuration  3 was  then  chosen  for  additional  modification  since  It 
had  the  lowest  flutter  speed  of  all  the  cases  possessing  strong  wing 
Interference  effects  (small  values  of  X). 

The  modifications  for  the  60  degree  sweep  case  Included  pitch 
and  roll  spring  stiffness  changes  and  the  addition  of  lumped  mass  to 
nodes  6 and  7 at  the  tip  trailing  edge  of  the  main  wing  (see  Fig.  12). 
Criteria  for  the  addition  of  the  lumped  mass  was  obtained  from  Reference 
12,  which  shows  the  general  effects  of  concentrated  mass  locations  on 
the  flutter  speed  of  swept  wings.  As  much  as  20%  of  the  total  wing  mass 
(10%  In  each  of  the  outer  two  sections)  was  added  during  the  analysis. 

Additional  configurations  were  Investigated  In  order  to  study 
the  effect  of  wing  sweep  on  the  flutter  speeds  of  the  modified  and 
unmodified  model.  The  configurations  that  were  evaluated  and  the  results 
obtained  are  given  In  Table  XV.  In  all  cases,  the  flutter  speeds  and 
frequencies  were  determined  by  the  k-method  of  analysis  using  the  zero 
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TABLE  XV 

l 

ANALYTICAL  FLUTTER  RESULTS 

No. 

Sweep 

uv 

— X 

“e  (ft) 

Modification 

VF 

ft/ sec) 

fF 

(Hz) 

1 

60° 

.8  .25 

None 

137 

6.8 

2 

60° 

.8  .90 

None 

140 

6.4 

3 

60° 

.6  .25 

None 

136 

(162)+ 

7.3 

(7.5) 

4 

60° 

.6  .90 

None 

136 

6.9 

5 

60° 

.4  .25 

None 

140 

6.0 

6 

60° 

.4  .90 

None 

142 

7.4 

7 

60° 

* .25 

Ml, (a.) 

124 

6.6 

8 

9 

60° 

* .25 

M2, (a.) 

121 

6.2 

- 

60° 

* .25 

R,(b. ) 

139 

7.0 

I 10 

11 

60° 

* .25 

P.(c.) 

119 

7.0 

60® 

.6  .25 

PM,(d.) 

107 

5.7 

12 

45® 

* .99 

None 

155 

7.2 

13 

45® 

# .99 

PM,(d.) 

107 

4.9 

14 

25® 

* 1.65 

None 

186 

8.3 

15 

25® 

# 1.65 

PM,(d.) 

150 

4.9 

X * Wing-Tall 

Horizontal  Separation 

* 

Same  fuselage  stiffness  as  configurations  No.  3 and  4 

* 

Same  fuselage  stiffness  as  configuration  No.  11 

+ 

Antisymmetric  aerodynamics  Instead  of  symmetric 

(a.) 

1.  10*  of  wing  mass  added  at  tip  trailing  edge 

2.  20*  of  wing  mass  added  at  tip  trailing  edge 

(b.) 

Roll  spring  flexibility  Increased  by  50* 

(c.) 

Pitch  spring  flexibility  Increased  by  50* 

<d.) 

50*  Increase  In  pitch  spring  flexibility  and  20*  of  wing  mass 
added  at  tip  trailing  edge 

1 

a 

1 

28 


value  of  structural  damping  as  the  Instability  point.  Computer  plots 
of  damping  versus  velocity  are  shown  In  Appendix  D for  the  configurations 
listed  In  Table  XV.  Also,  plots  of  the  flutter  modes  for  the  25,  45, 
and  60  degree  cases  are  given  In  Appendix  E for  both  the  standard  and 
modified  models.  It  should  be  pointed  out  that  these  plots  can  display 
the  character  of  the  flutter  eigenvector  at  only  one  Instant  of  time 
(t  * 0)  and  may  not  give  a complete  view  of  the  flutter  motion  and  pos- 
sible node  line  movement.  A true  picture  of  the  entire  flutter  mode 
can  be  obtained  by  displaying  the  real  part  of  the  complex  response  of 
each  structural  point  throughout  one  entire  cycle  of  motion. 

The  subsonic,  uncontrolled,  flutter  analyses  for  the  AFFDL 

model  were  carried  out  by  formulating  the  flutter  eigenvalue  problem  In 

12  14 

terms  of  the  natural  vibration  modes.  * The  current  study  employed, 

as  the  generalized  coordinates,  the  first  six  natural  modes,  which  were 

determined  by  the  methods  of  Section  II  and  Appendix  A.  The  flutter 

19 

equation  could  then  be  cast  In  the  following  form 

[K]-1  [[m]  + g(p)  [Q]]  (q>  = X (q>  (3.1) 

where 

[K]  Is  the  generalized  stiffness  matrix  (Ref.  Appendix  A) 
[m]  Is  the  generalized  mass  matrix  (Ref.  Appendix  A) 

[Q]  Is  the  complex  generalized  aerodynamic  coefficient 
matrix 

X Is  a complex  eigenvalue,  1 


r 


; 
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g(p)  is 


c*  Sp 


Cf  is  the  effective  root  chord 

S is  the  area  of  the  wing  semi  span 

p is  the  air  density 

{q}  are  the  normalized  complex  eigenvectors 


The  unsteady  generalized  aerodynamic  force  coefficients  used 

In  Eq.  (3.1)  were  calculated  by  the  doublet-lattice  method  of  Reference 

20,  which  was  developed  for  the  analysis  of  planar  Interfering  wing/ 

horizontal -tail  configurations.  These  aerodynamic  coefficients  were 

21 

then  modified  as  shown  by  Frederick  to  make  them  compatible  with  the 
form  of  Eq.  (3.1).  Before  computing  the  above  aerodynamic  force  terms, 
the  structural  modes  were  Interpolated  to  the  aerodynamic  grid  with  a 
local  least  squares  fitting  technique  developed  by  Cwach.1  The  aerody- 
namic modeling  for  the  25,  45,  and  60  degree  sweep  cases  Is  shown  In 
Figure  15. 

All  of  the  flutter  analyses  in  the  present  study  were  evalu- 
ated for  subsonic  flow  at  M s 0.4.  This  Mach  number  was  chosen  to 
encompass  all  possible  wind  tunnel  test  regimes  while  remaining  In  the 
region  of  fully  Incompressible  flow.  It  Is  shown  for  similar  analyses 
In  Reference  9 that  flutter  speed  variations  with  Mach  number  are 
negligible  below  M ■ 0.5. 

An  additional  consideration  In  the  flutter  analysis  Is  that 
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the  half-span  AFFDL  model  Is  suspended  such  that  It  simulates  anti- 
symmetric vibration  modes  and  mirror  Image  aerodynamics.  It  Is 
assumed  that  the  wind  tunnel  wall  represents  a reflecting  plane  through 
which  no  flow  penetrates.  This  requires  that  the  generalized  aerodynamic 
forces  be  calculated  assuming  a symmetric  loading  since  the  image 
system  is  performing  a symmetric  motion  to  satisfy  the  flow  boundary 
condition.  Due  to  the  antisymmetric  vibrational  characteristics  and 
the  symmetric  aerodynamics,  the  current  study  approximates  only  the 
antisymmetric  flutter  modes  under  wind  tunnel  test  conditions.  For 
comparison,  one  test  case  was  tried  with  antisymmetric  aerodynamics  and 
the  results  are  given  In  Table  XV  and  Appendix  D. 
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IV.  CONCLUDING  REMARKS 


The  results  of  this  study  have  demonstrated  very  good 
agreement  between  the  analytical  vibrational  modes  and  frequencies 
and  those  obtained  experimentally  for  the  60  degree  sweep  configuration. 
Additionally,  the  trends  of  the  analytical  vibration  studies  were  in 
general  agreement  with  the  results  given  In  References  1 and  9, 
Indicating  that  the  lumped  mass  dynamic  modeling  used  In  this  study 
was  quite  satisfactory.  All  differences  between  experimental  and 
theoretical  results  remained  within  5 percent. 

It  has  also  been  demonstrated  that  easily  made  modifications 
to  model  stiffness  and  mass  could  be  used  to  reduce  the  flutter  speeds 
and  frequencies  In  order  to  remain  within  the  model's  operating 
limitations  during  wind  tunnel  testing.  The  flutter  speed  for  the 
60  degree  configuration  was  reduced  as  much  as  21  percent  by  a combina- 
tion of  decreasing  the  pitch  spring  stiffness  and  adding  concentrated 
masses  near  the  tip  trailing  edge  of  the  wing. 

Additional  analysis  for  several  wing  sweep  configurations 
has  Indicated  that  flutter  speeds  of  the  model  can  be  Increased 
substantially  by  sweeping  the  wing  forward  to  25  degrees,  and  thus 
removing  It  from  close  proximity  to  the  tall.  Therefore,  the  emergency 
wing  sweep  mechanism,  which  rapidly  sweeps  the  wing  forward  on  command, 
should  be  effective  In  stopping  the  flutter  If  any  violent  Instability 
Is  encountered.  The  significant  variations  of  flutter  speed  with 
wing  sweep  were  attributed  to  changes  In  the  aerodynamic  Interaction 
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between  the  wing  and  tall  since  all  other  model  parameters  were 
constant  as  the  geometry  was  altered.  In  all  of  the  configurations 
analyzed,  the  mode  of  Instability  remained  the  same. 

When  using  the  results  of  this  study  to  determine  control 
laws  and  to  predict  the  wind  tunnel  behavior  of  the  actively  controlled 
model,  the  following  things  should  be  considered: 

1.  The  hydraulic  links  in  the  dynamic  model  were  considered  to 
be  rigid;  In  actuality,  they  possess  some  stiffness  charac- 
teristics which  were  not  determined  In  this  study. 

Additional  experimental  Investigation  of  the  model's 
dynamic  properties  with  a pressurized  control  system  is 
desirable. 

2.  The  aerodynamic  modeling  used  In  the  flutter  analyses 
assumes  smooth,  uniform  flow  over  the  wing  and  tall  and 
does  not  Include  the  effects  of  interference  between  the 
fuselage  housing  and  the  wing-tall  combination.  Studies 
of  the  flow  properties  around  the  model  for  a particular 
wind  tunnel  Installation  will  be  required  to  assess  the 
validity  of  these  assumptions  and  additional  aerodynamic 
fairings  may  be  required  to  produce  desirable  flow  qualities. 

3.  Previous  experience  has  shown  the  doublet-lattice  aerodynamic 
code  used  In  this  study  to  be  conservative  In  predicting 
flutter  speeds.  Also,  since  the  zero  value  of  structural 


damping  was  used  as  the  instability  point,  some  of  the 


Indicated  flutter  speeds  may  Increase  by  as  much  as 
5 ft/sec  when  the  true  structural  damping  is  used. 


I 
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Figure  1.  AFFDL  Model  Modified  for  Active  Stabilator  Control 


TUNNEL  CEILING 

TMTFPrnMWPrTnp  _ K A fiisfi  agf  torsion  spring 


Figure  2.  Schematic  of  Actively  Controlled  AFFDL  Wing-Tail  Flutter  Model 
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Figure  3.  Installation  of  Modified  AFFDL  Wing-Tall  Flutter  Model 
In  7 Ft.  by  10  Ft.  Subsonic  Wind  Tunnel 
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Figure  4.  Wing  Construction  of  AFFDL  Model 


Figure  8.  Schematic  of  Spar  Bending  Flexibility  Measurement 
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Figure  10.  Measurement  of  Spar  Bending  Flexibility 


Figure  11.  Measurement  of  Spar  Torsional  Flexibility 


Figure  12.  AFFDL  60°  Dynamic  Modeling 


SHAKER 


Figure  13.  Schematic  Diagram  of  Experimental  Modeshape  Measurement 
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TORSION  BAR  LENGTH  (IN) 

Figure  16.  Uncoupled  Fuselage  Torsion  Frequency  vs.  Torsion  Bar  Length, 
25-Degree  Sweep  Angle 


APPENDIX  A 


P 

COMPUTATION  OF  NATURAL  MODES  AND  FREQUENCIES 
FROM  EXPERIMENTAL  DATA 

Several  methods  are  available  for  calculating  the  generalized 
stiffness,  Krs,  and  generalized  mass,  mrs,  terms  appearing  in  the 
flutter  eigenvalue  problem,  Eq.  (3.1).  The  following  description  outlines 
the  method  employed  in  this  study  to  calculate  the  vibrational  modes  and 
frequencies  used  in  the  evaluation  of  these  terms.  Whenever  possible, 
experimental  measurements  of  mass,  stiffness  and  inertia  were  utilized 
In  the  analysis  so  that  the  calculated  natural  vibration  modes  and 
frequencies  would  closely  match  those  of  the  actual  model.  The  following 
analysis  procedure  was  Incorporated  into  a computer  code  which  would 

generate  the  lumped  mass  model  from  experimental  data  and  then  calculate 

♦ 

the  required  modes,  frequencies,  and  generalized  masses  for  a fully  flexible 
model.  The  computer  code  was  designed  such  that  variations  In  model  mass 
or  stiffness  could  be  easily  analyzed.  In  addition,  the  program  was 
written  so  that  It  could  be  extended  to  the  analysis  of  composite  wing 
and  tall  spar  configurations.  All  computations  utilized  the  classical 
vibration  and  modeling  techniques  employed  In  structural  dynamics  and 
aeroelasticity. 14,15,16 
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A.l  Lumped  Mass  Modeling  from  Experimental  Inertia 
and  Static  Unbalance 

Each  wing  and  tail  section  was  replaced  by  an  equivalent  dumb- 
bell unit  which  contained  equal  concentrated  mass  at  each  end  con- 

Q 

nected  by  a rigid,  massless  bar  (see  Fig.  12).  The  bar  was  positioned 
perpendicular  to  the  spar  center  line  at  the  spanwise  center  of  gravity 
location  for  the  section.  The  bar  lengths  forward  and  aft  of  the  spar 
center  line  were  then  calculated  for  each  section  such  that  section 
Inertia  and  static  unbalance  were  preserved.  This  calculation  is  shown 
in  Eq.  (A.l): 


where 

S Is  the  section  static  mass  unbalance  (slug-ft) 

I$  Is  the  section  Inertia  about  the  spar  line  (slug-ft2) 
M Is  the  total  section  mass  (slugs). 

The  wing  carry-through  assembly  was  also  modeled  with  a dumbbell  unit 
since  It  was  heavy  and  quite  flexible. 

Two  additional  concentrated  mass  points,  17  and  28,  account  for 
the  experimental  roll  inertias  of  the  forward  and  aft  fuselage  sections, 


respectively  (see  Fig.  12).  The  locations  of  these  nodes  were  chosen  to 
give  a useful  structural  deflection  point  and  their  mass  values  were 
then  calculated  to  give  the  correct  Inertia  about  the  roll  axis. 

A. 2 Assembly  of  the  Model  Flexibility  Matrix 

The  structural  flexibility  matrix  used  to  compute  the  natural 
vibration  modes  Is  assembled  from  the  experimentally  measured  uncoupled 
bending  and  torsional  influence  coefficients.  The  assembly  procedure 
computes  the  deflection  at  a lumped  mass  point  on  the  structure,  due 
to  a unit  load  at  any  other  point,  by  summing  the  contributions  due  to 
spar  bending,  spar  torsion,  fuselage  torsion,  roll  spring  flexure,  and 
pitch  spring  flexure.  The  moment  arm  lengths  about  the  different  tor- 
sional flexure  axes  are  computed  as  a function  of  spar  sweep  angle.  In 
this  manner,  the  uncoupled  flexibility  influence,  coefficients  of  the 
model's  elastic  components  can  be  easily  measured  experimentally  and  then 
assembled  to  obtain  the  entire  system  flexibility  for  any  geometry. 

It  Is  Implied  In  the  assembly  of  a structural  flexibility 

12 

matrix  that  the  structure  Is  constrained  at  some  point  to  prevent  the 
deflections  from  being  unbounded,  and  In  general,  the  constraint  can  be 
Imposed  at  an  arbitrary  point  on  the  structure.^6  The  model  In  this 
study  Is  constrained  at  mass  point  28  on  the  tall  spar  rigid  carry- 
through,  allowing  the  fuselage  torsional  degree  of  freedom  to  provide 
for  relative  motion  between  wing  and  tall. 
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A. 3 Solution  for  Vibration  Modes  of  a Structure  with 
Rigid  Body  Degrees  of  Freedom 

Using  the  flexibility  formulation,  the  free  vibration  inodes  of 
a constrained  structure  can  be  computed  from 

[K]"1  [M]  {h>  = J-  {h>  (A. 2) 

" wr 

where 

[M]  is  the  diagonal  lumped  mass  matrix 

[K]"1  is  the  flexibility  matrix  or  Inverted  stiffness 
matrix 

uy  is  the  natural  frequency  of  the  rth  mode 

(h>r  is  the  r*h  mode  shape  or  eigenvector. 

Since  the  AFFDL  model  Includes  a rigid  body  roll  degree  of 
freedom,  the  mass  matrix  In  Eq.  (A. 2)  must  be  modified  as  shown  in 
References  17  and  18  In  order  to  obtain  the  f ree-to-rol 1 , coupled,  anti- 
symmetric modes  of  the  entire  model.  This  modification  Involves  pre- 
and  post-  multiplying  the  lumped  mass  matrix  by  a constraining  matrix, 
[C],  which  preserves  angular  momentum  about  the  model  roll  axis.  The 
constrained  mass  matrix  Is  defined  as 


["']  • [C]T  [H]  [C]  . 


(A. 3) 
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For  the  case  of  rigid  body  roll,  the  constraint  matrix  is  defined  by: 


where 


[C]  « 

[I] 

*R 

(y> 


[I]  - f-  rn  (y}T  CM] 

lR 


(A. 4) 


is  the  identity  matrix, 

Is  the  total  mass  moment  of  Inertia  of  the 
structure  about  Its  rigid  body  rotation  axis, 
is  a vector  of  displacements  of  the  structural 
mass  points  due  to  a unit  rigid-body  motion. 
(Assuming  linearized  theory  and  using  a unit 
rotation,  this  vector  Is  the  vector  of  dis- 
tances of  the  mass  points  from  the  rigid  body 
roll  axis.) 


.The  constrained  mass  matrix,  [M* ] , Is  now  used  with  the  con- 
strained flexibility  matrix,  described  In  Section  A. 2,  to  solve  for  the 
relative  motion  vibrational  inodes.  (The  absolute  motion  of  the  flexible 
structure  In  any  mode  can  be  obtained  by  adding  the  proper  rigid  body 
component  to  the  relative  motion  as  shown  In  Reference  18.)  Substituting 
[M* ] for  [H]  In  Eq.  (A. 2)  gives  the  following  eigenvalue  problem  which 
can  be  solved  for  the  relative  motion  flexible  modes: 


ST*  [C]T  [M]  [C]  <hR> 


1 


(1), 


* <vr 


(A. 5) 
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After  obtaining  the  relative  flexible  modes,  (hR),  they  are  transformed 
back  to  the  absolute  (free-free)  modes  by 


I 


I 


{hA>r  « [C]  (hR>r  (A. 6) 

as  shown  In  Reference  17. 

The  solution  of  Eq.  (A. 5)  for  the  AFFDL  model  Involved  only 
27  degrees  of  freedom  since  the  28th  point  was  constrained.  The  [C] 
matrix  was  computed  with  dimensions  of  28  x 28,  but  the  28th  column  Is 
dropped  during  the  calculation  of  Eq.  (A. 3)  in  order  to  reduce  out  the 
mass  at  the  constrained  point  (I.E.,  [C]  Is  used  as  a 28  x 27  matrix). 
The  entire  28  x 28  [C]  matrix  Is  used  In  Eq.  (A. 6)  to  transform  from 
the  relative  to  the  absolute  modeshapes. 

The  rigid  body  mode,  which  Is  a gravitational  pendulum  mode, 
can  now  be  calculated  from  one  of  the  relative  flexible  modes  by  the 
method  given  In  Reference  18.  In  matrix  form,  this  calculation  can  be 
written  as 


where 


(hA> 


rigid 


(y)  • 4>r 


(A. 7) 


<(>R  Is  a scalar  representing  the  value  of  the  rigid 
body  coordinate  (In  this  case  a rotation  angle). 

The  rigid  body  component,  tR,  can  be  expressed  as 


53 


♦R  “ - r {y>T  [M]  {hR}  (A. 8) 

R 

where 

{hR}  is  any  one  of  the  previously  determined 
relative  flexible  modes 

After  combining  Eq.  (A. 7)  with  Eq.  (A. 8),  the  rigid  pendulum  mode  is 
calculated  from 

{hA>  - - r-  (y)  {y}T  [M]  {h  > (A. 9) 

M rigid  *R  r 

Eq.  (A. 9)  Is  recognized  as  the  second  expression  of  the  constraining 

equation,  Eq.  (A. 4),  multiplied  by  a flexible  modeshape.  The  first 

flexible  mode  Is  used  In  this  study  since  It  Is  determined  with  the 

greatest  computational  accuracy. 

The  complete  set  of  free-free  vibration  modes  Is  now  normalized 

so  that  each  modeshape  vector  has  unity  magnitude.  It  should  also  be 

mentioned  at  this  point  that  the  solution  procedure  Insures  that  the 

flexible  modes  are  orthogonal  to  each  other  with  respect  to  the  lumped 

mass  matrix,  [M].  Additional  orthogonality  conditions  between  the  rigid 

11  u 

body  mode  and  the  flexible  modes  are  also  satisfied.  * 

Now,  the  generalized  mass  can  be  computed  from  the  absolute 
modes  of  Eqs.  (A. 6)  and  (A. 9)  by 


I 


" T!1*  -V 
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where 


"Vs  ■ {hA,Tr  W fhA>s 


{h«}  Is  the  vector  of  displacements  of  the  sth 
M s 

absolute  mode 


{h.}^  Is  the  transpose  of  the  vector  of  displace- 
M r 

ments  of  the  rth  absolute  mode. 

Because  of  the  orthogonality  condition  of  the  natural  modes  with  respect 
to  the  lumped  mass  matrix 

m^  ■ 0 for  r i s . 
rs 

The  generalized  stiffness,  Kr$,  Is  then  determined  from 


AFFDL  60°  Modeshapes  , — = 0.4 


Mode  2 a)  = 4.3  Hz 


Mode  3 w = 9.8  Hz 


APPENDIX  C 

EXPERIMENTAL  VIBRATION  MODES 


APPENDIX  D 


PLOTS  OF  DAMPING  VS.  VELOCITY  AND 
FREQUENCY  VS.  VELOCITY+ 


t For  abbreviations  and  symbols,  reference  Table  XV. 
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